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We present a new model for generalized parton distributions (GPDs), based on the aligned jet 
model, which successfully describes the deeply virtual Compton scattering (DVCS) data from HI, 
ZEUS, HERMES and CLAS. We also present an easily implementable and flexible algorithm for 
their construction. This new model is necessary since the most widely used models for GPDs, which 
are based on factorized double distributions, cannot, in their current form, describe the DVCS data 
when employed in a full QCD analysis. We demonstrate explicitly the reason for the shortcoming 
in the data description. We also highlight several non-perturbative input parameters which could 
be used to tune the GPDs, and the i-dependence, to the DVCS data using a fitting procedure. 



PACS numbers: ll.lO.Hi, ll.30.Ly, 12.38.Bx 



I. INTRODUCTION 

Generalized parton distributions (GPDs) have been 
studied extensively in recent years [Q, ||, ||, ||, ^, |^, 
^, H, |lo[ |ll|, |l^. This interest was spurred by the re- 
ahzation that these distributions are not only the basic, 
non-perturbative ingredients in hard, exclusive processes 
such as deeply virtual Compton scattering (DVCS), or 
exclusive vector meson production, but that they are gen- 
eralizations of the well known parton distribution func- 
tions (PDFs) from inclusive reactions. GPDs incorporate 
both a partonic and distributional amplitude behavior 
and hence contain more information about the hadronic 
degrees of freedom than PDFs. In fact, GPDs are true 
two-parton correlation functions, allowing access to the 
highly non-trivial parton correlations inside hadrons [ p3| . 

GPDs can be broadly characterized by the following 
features: 

• They depend on two momentum fraction variables, 
a partonic variable defined with respect either to 
the incoming or the average of the incoming and 
outgoing proton momentum and the skewedness 
(which is the difference between the momentum 
fractions of two adjacent partons in the parton lad- 
der). 

• For fixed skewedness, they are continuous functions 
of the dependent variable and span two distinct re- 
gions, the DGLAP region and the ERBL region, 
in which their evolution in scale obeys generalized 
versions of the DGLAP and ERBL evolution equa- 
tions, respectively, and in which their behavior is 
qualitatively different. 



• They are even functions of the skewedness variable 
and the singlet, non-singlet and gluon distributions 
are either symmetric or anti-symmetric about the 
center point of the ERBL region (the symmetry 
obeyed depends on the precise definitions used). 

• The Lorentz structure of their definitions implies 
a polynomiality condition |^, |ri|, |TJ]: their 
{N — l)-th moments are polynomials in square of 
the skewedness of degree no greater than N/2. 

• They reduce to the ordinary PDFs in the limit of 
zero skewedness (the 'forward limit'). 

All of the above features have to be preserved under evo- 
lution in scale. 

Any suggested model of GPDs should adhere to these 
mathematical features. In jl^ such a model, based on 
double distributions (DDs), was suggested for the GPD 
input distributions (see also 0). In Q it was real- 
ized that an additional term, the so-called D-term, was 
required in the ERBL region for the unpolarized quark 
singlet and gluon distributions in order to satisfy poly- 
nomiality for even A''. The use of factorized DDs 
augmented with a D-term has become a popular phe- 
nomenological model. Unfortunately, when this type of 
model for input GPDs was used in its current form to cal- 
culate deeply virtual Compton scattering at both leading 
(LO) and next-to- leading order (NLO), the results were 
not in agreement with the HI data ||l8| on the DVCS pho- 
ton level cross section, a{'y*p jp), and the HERMES 
and CLAS data jl^ on the DVCS single spin asymmetry 
or charge asymmetry [^0|, p^ . 

Another popular model for input GPDs, inspired by 
the aligned jet model (AJM) |Q and its QCD exten- 
is based on the observation that at a scale 
2 GeV^ and a wide range of x^j , soft physics 
gives the dominant contribution to the parton densities. 
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As a result the effect of skewedness at small Xbj should be 
rather small and hence at the input scale, it is a good ap- 
proximation to set the GPDs equal to the forward PDFs 
at the same parton fraction, X, defined with respect to 
the incoming proton |^ (for any skewedness). This has 
the advantage that it automatically satisfies the require- 
ments of polynomiality for the first two moments, how- 
ever one encounters infinities in the quark singlet GPD 
in the middle of the ERBL region. 

Another 'forward model', which may be considered to 
be an extreme case of a DD model, was adopted in Q 
where one assumes that the GPD is equal to the forward 
PDF at the same parton momentum fraction, v, with re- 
spect to the average of the incoming and outgoing proton 
momentum P = [p + p')/2 which implicitly contains the 
skewedness. This translates to an X which is shifted to 
lower values by an amount controlled by the skewedness. 
This ansatz works fine for the DGLAP region. Unfortu- 
nately in the ERBL region it also involves sampling the 
forward PDFs right the way down to zero in momentum 
fraction where they have not yet been measured (this is 
especially problematic for singular quark distributions). 
In this paper we construct an alternative, finite 'forward 
model' for the input GPDs, using the forward input PDFs 
in the DGLAP region and imposing a simple form in the 
ERBL region that has the correct symmetries and ensures 
pol ynomiality is respected in the first two moments (see 
for alternative ways of dealing with this problem). 
As we will demonstrate this forward model reproduces 
the available data on DVCS reasonably well. 

Throughout this paper we will use the off-diagonal rep- 
resentation of GPDs, J-^{X, C), defined by Golec-Biernat 
and Martin [ p^ and used, for example, in the numeri- 
cal solution of the GPD evolution equations in ||2^ (see 
pO , 26, 2^ for other approaches to numerical evolution). 
They depend on the momentum fraction X G [0, 1] of 
the incoming proton's momentum, p, and the skewedness 
variable ( = A+/p+ (so that C = for DVCS). For the 
quark case, the relationship of the quark and anti-quark 
distributions, (X , C) , J^H^ X) , to Ji's GPD H'J{v,^) 
is shown in Fig. |l| with ^ = A+/2P+. More explicitly, 
for V e [~^, 1]: 



and for v G [—1, £\ 



l-C/2 



l-C/2 



(1) 



(2) 



The two distinct transformations between v and X for 
the quark and anti-quark cases are shown explicitly on 
the left hand side of Eqs.(|l|, H) . There are two distinct 
regions: the DGLAP region, X > C, {\v\ > ^), and the 
ERBL region, X < C, {\v\ < C) . In the ERBL region, due 
to the fermion symmetry, J-'^ and T'^ are not indepen- 
dent. In fact JP«(X, C) = -T^^iC -X,C), which leads to 
an anti-symmetry of the unpolarized quark singlet dis- 
tributions (summed over fiavor a), J"^ = EaT^^" + T"^'", 
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FIG. 1: The relationship between T'^{X, (), T'^{X, C) and Ji's 
function H''{v,^) with v G [-1, 1] and X G [0, 1]. 
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FIG. 2: Symmetric double distributions (left), indicating mo- 
mentum fractions of the outgoing and returning partons, and 
(right) their physical domain. 



about the point C/2 (the non-singlet and the gluon, JTS, 
which is built from vHj^{v,S,), are symmetric about this 
point). J^'^^'S are constructed analogously to from 

The rest of the paper is organized as follows. Section || 
contains a detailed explanation of why DD-based models 
in their current form cannot describe the data. In Sec- 
tion III we construct our alternative forward model for 

~ M and 



input GPDs, which is motivated by the AJM 
describes the data well. In Section IV we propose a phe- 
nomenological model for the slope of the i-dependence in 
which the slope parameter is allowed to change with pho- 
ton virtuality, — —Q^. The model improves the the- 
oretical description of the Q^-dependence of the HERA 
data, relative to using a constant slope. Finally we sum- 
marize our findings in Section |^. 



II. THE PROBLEM WITH DOUBLE 
DISTRIBUTION MODELS 

In this section we discuss factorized double distribution 
(DD) based models for GPDs and explain in detail why 
the sampling of the forward PDF at extremely small x 
in constructing the GPD leads to a problem in the quark 
singlet GPD. 
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Symmetric DDs, FDD{x,y,t,Q^), were introduced in 
with plus momentum fractions, x, y, of the outgoing 
and returning partons defined as shown in the left hand 
plot of Fig. y. They exist on the diamond-shaped domain 
shown to the right of Fig. ||. The outgoing parton lines of 
course only have a single plus momentum relative to any 
particular external momenta, so the GPDs are related to 
these DDs via a reduction integral, involving S{v — x—^y), 
along the off- vertical lines in the diamond (the dotted line 
corresponds to w = ^): 

H{v,0^ dx dy'6{x' + S,y' - v)Fdd{x' ,y'). 

J -I J -l+\x'\ 

(3) 

In Ijlfl, 17| a model for F£i£i{x,y,t, ^'^), at the input 



the effects of skewed evolution are much stronger in the 
region of X ~ compared to forward evolution ^ |5) . 
Note that this region strongly influences the cross section 
at small Xh-\ and some asymmetries both at small and 
large Xhj p9| , if they are dominated by the imaginary 
part of the amplitude which is proportional to T^((,,C) 
at LO (see for example 21 as well as pO, Bl^, p2]). 



One may wonder whether one can come closer to the 
data by choosing a very low input scale and valence-like 
partons as in the GRV scenario generating the rise 
of the parton distributions entirely through evolution. It 
turns out that choosing the canonical value of 6g = 1 and 
GRV98 input distributions (Qo = 0.51,0.63 GeV in LO 
and NLO, respectively) the curves still overshoot the data 
considerably. Indeed, even if one tries to minimize the 
effect of the enhancement due to skewedness at the input 
Ql , was introduced\n which' the functional ^cale, by choosing a large value of 6, = 100, evolution still 



form is factorized in x and y and in the i-dependence, an 
assumption borne out of convenience rather than physical 
necessity: 



F\ 



DD(x,y,ix\t) = n'{x,y) f (x,/)r*(t). 



(4) 



Here and are the standard PDF and form factor 
for the parton distribution of general type i. Since t- 
dependence is assumed to factorize and thus of no im- 
portance to the following, we will suppress it from now 
on. The profile functions, 7r*(a;,y), are asymptotic shape 
functions [0 for quarks and gluons of the general form 



Tr{x,y) 



r(26 + 2) [il-\x\)^-y 



2lb 



22b+ir2(5+i) (i_|a;|)2fc+i 
and normalized such that 



(5) 



drives the prediction above the data by at least a factor 
of ^ 4. Since, for example, daovcs — |-^'^(CiC)p7 the 
large enhancement of the quark singlet GPD &i X — C, 
is the root of the problem in this model. As we will 
explain below the origin of the enhancement stems from 
sampling singular forward sea distributions at extremely 
small X in the DD-based model. To understand the last 
statement, one has to first establish the regions in x in 
which the PDFs are sampled in the double distribution 
model, particularly at small x. 

First, having defined the model for the factorized dou- 
ble distribution in Eq. (^) one may then perform the y'- 
integration in Eq. (||) using the delta function. This then 
modifies the limits on the x'-integration according to the 
region concerned: for the DGLAP region X > C, {v > S^) 
one has for the quark GPD 



l + \x\ 



dy n{x,y) = 1 , 



(6) 



Note that tt is an even function of both of its arguments 
p8[ . The power b controls the size of the skewing ef- 
fects in the input GPD. Usually 6 = 1 is chosen for the 
quarks, corresponding to maximum skewedness, whereas 
6 = 2 is chosen for the gluons. In the limit 6 — > (X) there 
is no external skewedness effect, however, since in this 
limit, H{v,£^) — q{v), this translates into internal depen- 
dence on the skewedness J?^(X, (^) = q /(l~C/2)- 

Note that a consequence of the above model is a ratio of 
GPD/PDF for the quark singlet in the DGLAP region at 
X = ( = Xbj , which is substantially larger than 1 |g5[ |2^ 
for all experimentally relevant values of xi,j . 

Concerning the description of the data using the above 
model, it was shown in pl} | that maximal skewing (b = 1) 
at 'conventional' input scales (Qo = 1,2 GeV) overshoots 
the HI data by a factor of 6 — 10. It was also demon- 
strated that one can describe, in LO only, the HI data 
without including skewedness effects at the input scale, 
Qo = 2 GeV, if one neglects evolution |2|]. This simpli- 
fication, however, is not warranted since we know that 



T^'^{XX) = 

^ ^ dx'^" [x', -^{X- x') + x'-l^ q%x') . (7) 

For the anti-quark GPD in the DGLAP region X > 
C (v < — ^) one has 

T^''^{X,0 = 

^ J dx'n" (x', ~^{X + x') + x' + lj it{\x'\) , 

(8) 

changing variables from x — » —x and exploiting the fact 
that the profile functions are even in both arguments one 
arrives at 



T^-'^iXX) 



I ^ dx'n" (^x', |(X ~ x') +x'-l^ r{\x'\) , (9) 
so that the singlet and non-singlet quark distributions 
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are given by 

a 



X 



where y{x') = 2{X - x')/C + x' ~l. 

In the ERBL region, X < {\v\ < ^) integration over 
y leads to: 



f dx'iT'' (^x',^{X ~x')+x' -l^q'^ix') 



ix-c 
T^'^iXX)^ 



dx'TT'J (^x',-^{X + x') + x' + 1^ q'^ilx'l) 
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Jo 
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dx'TT'J 


J-x 




jain, 


using X 











dx'n" ( x', --{X + x') + x' + 1] q''{x') - 



(11) 
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dx'n'^ (x , y{—x')) <t{x') 
C-x 

dx'ir'i {x',y{-x'))q''{x') 



dx'Ti" {x\y{x'))q''{x') 



(12) 



Hence, for the singlet and non-singlet combinations one 
has 



Y.\[[ dx'n^ {x',y{x')) [q^{x') + g-(x')] 

C-x 



dx'n^ix^yi-x')) [q-{x') + q-{x')] 



'0 



dx'TT^x^yix')) [q-ix')-r{x')] 



K-x 



+ / dx'nHx\y{-x')) [q'^{x')-q'^{x') 
Jo 



(13) 



These expressions clearly satisfy the correct symmetries 
properties, i.e., i^^(C - ^, C) = -^"^(^,0, F^^-''{C - 
XX) = F^^'"(X,C) (n.b. y{x') y{~x') when X -> 
C — X ). Analogously for the gluon one obtains 

J^^iX, 0^7 r '^^'^^ (^'' ' (14) 

for the DGLAP region, and 

^^iXX) = ^\ r dx'n^x\yix')) x' g{x') 
^ '-Jo 

K-X 

+ / dx'-K^ {x',y{-x')) x'g{x') , (15) 
for the ERBL region (which is symmetric under X 

C-X). 

Inspection of the integration limits in Eqs. (|l^, [l^ ) 
highlights the main problem. In the limit A" — > ^, as a 
result of the lower limits of the integrals the forward PDF 
is sampled closer and closer to x' = 0, where it has not 
yet been measured. This will be irrelevant providing the 
integrand is sufficiently non-singular in x' in this region 
i.e. it can happen if the profile functions, tt', provide a 
strong suppression of this region, or if the PDFs them- 
selves are sufficiently non-singular. However, we know 
that phenomenological quark and sometimes even gluon 
input distributions are singular in the small x region. 
In the quark case this problem is made worse by the 
fact that we sample the number distribution q{x') rather 
than the momentum distribution x' q{x') (so that a non- 
singular momentum distribution xq{x) oc a;" for a G [0, 1] 
will give a singular number distribution q{x) oc a;""^). 

It turns out that for realistic quark distributions the re- 
gion close is x' = is very significantly sampled for small 
C = Xbj- This leads to two serious problems. Firstly, 
the forward distributions are unknown here so one must 
extrapolate the 'known' analytic forms downwards in x' . 
Secondly, and much more importantly, it leads to a very 
significant enhancement of the quark singlet GPDs rela- 
tive to the PDFs for AT w C, i.e., the region most relevant 
for DVCS. Though of paramount importance for DVCS, 
this region is but a small region of phase space where the 
current factorized DD models fail. 

We illustrate this using a series of three figures relating 
to the formation of the quark singlet GPD in the DGLAP 
region close to X = (. Fig. shows the integrand, 
I{x') = 7r"(a;', y) u^{x'), of Eq.(pX|) for the up quark sin- 
glet (multiplied by C) as a function of x' for two values 
of C = 0.1, 0.0001 and two values of AT - C = 0.1, 0.001. 
Clearly as X approaches ( the PDF is sampled at pro- 
gressively smaller values of x' <^ C, where for small C. it is 
unknown. Fig. ^ shows the average value of x' sampled in 
this integral (divided by () as a function of ( for several 
values oi X — ^. For very small values of X — ^ the average 
value of x' settles down to about C/4, for small ^. Finally 
in Fig. H we show, for MRST input PDFs, the ratio of 
the quark singlet GPD to PDF at C = 0.0001, 0.1, for the 
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FIG. 3: The integrand of Eq.(p^, illustrating how the up 
singlet PDF is sampled in the DGLAP region close to the 
boundary of the ERBL region, to produce the up singlet GPD. 




FIG. 4: The average value of x' sampled in the DGLAP region 
in the double distribution model, for the up singlet GPD, close 
to the boundary with the ERBL region as a function of the 
skewedness. Several values of X ~ ( are shown. 



canonical value of the power, bq = l, inEq. (||). Note the 
large enhancement of the GPD at X « C, particularly for 
small C in the upper plot. We emphasize that this en- 
hancement, which leads to an overshoot of the DVCS 
data, is built in right at the start in the modeling of the 
quark singlet GPD at the input scale. One also sees that 
for the gluon, which uses x'g{x'), and bg = 2, the ratio 
remains close to unity. 

The most important enhancement effect in the valence 
region, C > 0.1, originates from the relative shift of the 
parton momentum fraction X to smaller values close to 
X — (although the enhancement from small (x ) is 
still significant). As we will show in the next section, 
the assumption that H{v,^) = (1 - C/2)T{X,C) = q{v) 
with V = {X — C/2)(l — C/2) gives a good description of 
the data at both small and large Xbj- As stated before. 
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FIG. 5: The ratio GPD to PDF at C = 0.0001 (upper plot) 
and = 0.1 (lower plot) for the quark singlet and gluon in 
the double distribution model, using MRSTOl distributions 
in LO and NLO, at the input scale Qo = 1 GeV. Note the 
large enhancement of the quark singlet close to X = (. 



this corresponds to a factorized DD-model with b — oo, 
i.e., with no external skewedness. However, in terms of a 
comparison of GPD to forward PDF, there is a residual 
effect of skewedness since one now has to compare q{v) 
with q{X). Since we are comparing number distributions 
which are more singular than momentum distributions, 
any shift in the momentum fraction to smaller values 
will lead to a quite a large enhancement of q{v) relative 
to q{X). For CTEQ6M for example the enhancement at 
X = 0.1 and ( — 0.1 is about 1.7 for the quark singlet, 
which increases further if more skewedness is added by 
decreasing b. 

However, as we will demonstrate in Section HI, the 



available data allows little room for further enhancement 
due to skewedness at the input scale since the LO result, 
at least, is already close to the upper bound of the ex- 
perimental errors. Therefore, only the extremal "6 — cxd" 
version of the current factorized DD-model can be used 
to describe the data. An obvious solution to this is to 
modify the quark singlet profile functions in Eq.(M) in 
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such a way as to suppress the region of very small x' . 
However, one must find a new functional form which 
achieves this without spoiling the known mathematical 
features of GPDs discussed above. Alternatively, one 
abandons the attempt to model a DD using a factorized 
form Tr{x,y)f{x), though appealing due to its simplicity 
but possibly too simplistic in its form. This remains an 
open problem and has to be addressed by those who wish 
to use the double distribution framework, factorized or 
not, to model GPDs. 
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III. THE FORWARD INPUT MODEL AND THE 
ALIGNED JET MODEL 



In this section we revisit the logic of setting the 
GPDs equal to the forward PDFs by proposing an al- 
ternative forward model to that suggested in ||ll|, with 
suitably-symmetrized input GPDs in the ERBL region 
constructed so as to satisfy the requirements of polyno- 
miality for the first two moments. 

In [pi DVCS was predicted to be measurable at DESY- 
HERX and, allowing for the freedom associated with 
choosing the slope parameter, B, the predictions success- 
fully describe both the HI data |18 and the recent ZEUS 



result 1 34, ^ on the photon- level DVCS cross section. 
This was achieved by modeling the imaginary part of 
the DVCS amplitude at the input scale using the aligned 
jet model (AJM) iQ. This was then compared to the 
imaginary part of the DIS amplitude, calculated within 
the same framework, which was found to be smaller by a 
factor of about two. The comparison enabled the normal- 
ization of the DVCS amplitude at the input scale to be 
set using F2 structure function data. The DVCS ampli- 
tude was then evolved to higher scales using LO skewed 
evolution in perturbative QCD. 

The basic relation between the DVCS and DIS ampli- 
tudes, using the AJM, is given by M 



R 



IniTr 



DVCS 



IniTr 



In 1 



DIS 



91 

Mi 



Ml 



1.5 - 2.5, 
(16) 



where for the AJM is typically 1-3 GeV^ and Mq is 
a hadronic scale which roughly corresponds to the low- 
est allowed, excited intermediate state in the s-channel. 
Therefore, M§ - 0.4 - 0.6 GeV^ or about m^. The 
AJM neglects the contribution of quarks with large trans- 
verse momenta in the quark loop attached to the pho- 
tons in the handbag diagram. Since the contribution of 
small transverse momenta is more symmetric than the 
one at large transverse momenta, the AJM may some- 
what overestimate the effect of skewedness at the input 
scale. Eq. ( p^ can be generalized to demonstrate how 
the forward limit ImTpyQg = ImTpjg is achieved, i.e., 
how the skewedness effect is reduced by giving the out- 
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FIG. 6: The ratio ii as a function of X = Q ^/Q^ for several 
value of and two values of Mq — 0.4 GeV^ (upper plot) 
and 0.8 GeV^ (lower plot). 



going photon a space- like virtuality, = —Q 



ImTr 



DIS 



1 







MS 

w 









(17) 



This procedure allows one to derive a very important 
relation between the relative momentum fractions of the 
outgoing and returning partons, X and —{X — C) of the 
quark singlet GPD, and the virtualities of the incoming 
and outgoing photons 



A = 



X Q- 



(18) 



We illustrate this point in Fig. ^ by plotting i? as a 
function of A for several values of and two values of 
Mq to demonstrate the relative insensitivity of R (within 
20 ~ 30%) to Mq. 

The plot shows that as X increases relative to C in 
the DGLAP region the ratio drops rapidly to its forward 
limit. For example, at A = 1/2, i.e., X = 2^, the curves 
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are very flat and there is only a modest enhancement of 
20 ~ 40%. One also encounters this behavior in the DD 
model if one investigates the ratio of the GPD to the PDF 
in the DGLAP region (see for example ||2^, ^). It would 
therefore be advantageous to be able to directly relate R 
to a ratio of GPD to PDF. Trusting that perturbative 
QCD is applicable at the AJM input scale one can, in LO 
at least where the coefficient function is trivial, directly 
translate the ratio in amplitudes for a particular A into 
a ratio of GPD to PDF: 



2.4 



R{\) = 



IniTr 



DVCS 



ImT^jg 

7S (V - 



(X^C/(1-A),C) 
{l~C/2)qS{X) 



I.e., 



(C/(l-A),C)-(l-C/2) i?(A) q^X). 



(19) 



(20) 



There are several comments in order at this point: A is 
now bounded from above through A < 1 — C- This im- 
plies that the relationship between the ratios in Eq. ( |20| ) 
is only strictly true for A = (i.e., for DVCS for which 
— Q). The case A ^ should be viewed as follows: for 
^ 2 GeV^, there is still the possibility of having more 
than one rung in the partonic ladder. Probing the upper- 
most rung with a virtuality reveals the distribution in 
momentum fractions, in this case X = C, i.e., X — C, = Q 
corresponding to A = 0. The next rung and its distri- 
bution in momentum fractions can be accessed b^ 'emit- 
ting' a photon with space-like virtuality (i.e., > 0). 
As Q ^ increases and one goes further down the ladder 
to where X ;2> one approaches the forward limit. If 
one keeps the interpretation of the s-channel cut as be- 
ing equal to the imaginary part of the 'scattering' ampli- 
tude for Q ^ ^ 0, which, in LO, is directly proportional 
to GPD/PDF at X 7^ C rather than at X = then 
the ratio i? is a direct measure of quark singlet GPD to 
the quark singlet PDF for X ^ C,. However, this logic is 
valid only at LO. At NLO, the situation radically changes 
since, first, the gluon directly enters into the amplitude, 
secondly, the convolution of the coefficient function with 
the GPD is no longer as trivial as in LO and thirdly, the 
value of Us at low is quite different in LO and NLO. 
Therefore such a simple relation as in Eq. (20) should 
and can no longer be valid. In order to keep the ansatz 
as simple as possible, we will only require that the model 
GPDs, at least in LO, produce a ratio, i?, which is in 
broad agreement with the R values obtained in the AJM 
from Eq. (|20|). 

If one chooses the forward model ansatz where the 
GPD equals the PDF at v in both LO and NLO, due 
to a lack of a better Ansatz, (see e.g. pj): 



q^{v) EE q^ 



X - 



C/2\ 



1 - c/2 ; 



(21) 



which corresponds to the b ^ oo limit of the DD model, 
one obtains a ratio of GPD to PDF at X = C of ~ 2.3 
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FIG. 7: The ratio GPD to PDF at C = 0.0001 for the quark 
singlet and gluon, using MRSTOl distributions in LO and 
NLO, at the input scale Qo ~ 1 GeV. This ratio is weakly 
dependent on for small i^. 



for the quark singlet, in agreement with the AJM pre- 
diction though for slightly different values of Qq (di- 
rectly compare the upper line in the upper plot of Fig. ^ 
with the quark singlet in Fig. ^ keeping in mind that 
A = 4=> X/C = 1 and A = 0.9 <^ X/C = 10). In con- 
sequence, our model ansatz corresponds to an AJM with 
maximal skewedness. 

For our forward model in the DGLAP region we, there- 
fore, choose for simplicity, the ansatz of Eq. ( pi] ) for the 
quark singlet, the non-singlet (i.e., the valence) and the 
gluon. The ratios of GPD to PDF at the input scale 
for MRSTOl quark singlet and gluon distributions |3^] 
at LO and NLO at Qo = 1 GeV are shown in Fig. |^. For 
a function which falls as x decreases, such as the valence 
quark at small x or the MRST gluon at LO, this ansatz 
leads to a suppression of the point X = relative to the 
forward case (see the lower, dotted line in Fig. 0). Note 
that in NLO the MRST gluon actually goes negative at 
small X, so a ratio of GPD/PDF > 1 close to X = C in 
this case leads to a suppression of the DVCS cross sec- 
tion from the gluon contribution, relative to using PDFs. 
However, the DVCS cross section is rather insensitive to 
the behavior of this ratio close to X = for the gluon, 
since it only enters in NLO, and is completely insensitive 
to it for the non-singlet quark case since this distribution 
only enters into the evolution. 

The above reasoning indicates that the physics of the 
AJM model provides a guide for modeling input GPDs in 
the DGLAP region. If one compares the NLO imaginary 
part of the DVCS amplitude from the above model to the 
NLO imaginary part of the DIS amplitude as extracted 
from a recent HI F2 fit @, we find for MRST2001, R ~ 
3.8 and for CTEQ6M, R ~ 2.7 for = 3 GeV^ and x = 
0.0005. the enhancement effect of evolution, the AJM 
result for R is basically reproduced for CTEQ6M but not 
for MRST2001 using eq. (Ell). Given the widely different 



8 



parameterizations at NLO, this seems acceptable to us, 
at present. As we will see below the enhancement effect 
generated through the shift is too strong for both LO and 
NLO at low values of near the input scale. 

The prescription in eq. docs not dictate what to do 
in the ERBL region, which does not have a forward ana- 
log. Naturally the GPDs should be continuous through 
the point X = C, and should have the correct symmetries 
around the midpoint of the ERBL region. They are also 
required to satisfy the requirements of polynomiality: 

Mm = j dvv^-^[H'i{v,^)+HS{v,0] 

N/2 

= ^C'"C2fc,iv. (22) 

k=0 

At this point we choose to model the ERBL region 
with these natural features in mind. We demand that 
the resultant GPDs reproduce the first moment Mi « 3 
and the second moment M2 « 1 (with the D-term set 
to zero to remove the quadratic piece in Eq. (Bil)) where 



Ml 
M2 



1 



dwiI«'^'^(w,C) =3, 



(23) 



(24) 



and C in Eq. (|24|) was computed in the chiral-quark- 
soliton model p8|| . This reasoning suggests the following 
simple analytical form for the ERBL region {X < (): 



jrf'^^(x,c) = J^^'^^(C) [1 
:f^xx)^J'HO 



C/2 



AS'^'^iOC^'^^iXX)] 
[l + A'^iOC'^XX)] 



(25) 



where the functions 



c^^^^ix,c) = l-^ 



c^ixx) 



15 

y 




(26) 



vanish at X = ^ to guarantee continuity of the GPDs. 
The A*(^) are then calculated for each ( by demanding 
that the first two moments of the GPDs are explicitly sat- 
isfied (remembering to include the D-term in the ERBL 
region which only provides the quadratic term in ^ in 
Eq. ([23|)). For the second moment what we do in prac- 
tice is to set the D-term to zero and demand that for 
each flavor the whole integral over the GPD is equal to 
the whole integral over the forward PDF for the input dis- 
tribution concerned (due to the inherent small errors on 
the PDFs, the sum of such integrals will be close to, but 
not precisely equal to, unity). Note that the modeling 
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of the ERBL region is unimportant at small ( since the 
unknown subtraction constant in the dispersion relation 
between the real part of the amplitude which formally 
depends on both the ERBL and the DGLAP region, and 
the imaginary part which formally depends only on the 
DGLAP region, is proportional to ( and therefore incon- 
sequential at small C.. 

It would be straightforward to extend this algorithm 
to satisfy polynomiality to arbitrary accuracy by writing 
the A*(C) explicitly as a polynomial in ( where the first 
few coefficients are set by the first two moments and the 
other coefficients are then randomly chosen since nothing 
is known about them. The above algorithm is extremely 
fiexible both in its implementation and adaption to ei- 
ther other forward PDFs or other functional forms in the 
ERBL region. Therefore it can be easily incorporated 
into a fitting procedure. 

In Fig. ^ we show the shape of the resulting input 
GPDs for two characteristic values of C = 0.001, 0.1. The 
upper plot in this figure explicitly shows the antisymme- 
try of the singlet GPD and the symmetry of the gluon 
GPD about the point X = C/2. 

The photon level cross section results from this model, 
using MRSTOl [|| and CTEQ6 ^ input distributions 
at LO and NLO, are compared in Figs. §|l^to the HI |l|] 
and ZEUS p5| data at their average kinematic points, 
respectively. In these curves we chose to use an x and 
Q^-independent slope parameter of B — 6.5 GeV~^, but 
realistically there is a 30 — 40% uncertainty associated 
with the value of this unknown parameter. The figures 
illustrate that within the framework of the forward input 
model for GPDs the DVCS cross section remains rather 
sensitive to the choice of input PDF and to the accu- 
racy with which the calculation is performed (i.e., LO or 
NLO). It should be noted that the description at lowest 
HI value of is bad. However, the enhancement we 
chose using eq. ( |2l] ) at Qq = 1 — 1.69 GeV^ corresponds 
more to the AJM at = 3 GeV^ (compare Figs. 
Hence it is not surprising that the description at low 
is not good suggesting that the shift in X in eq. ( |2l| ) 
should be less at lower values of Qq . When increasing 
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FIG. 8: The quark singlet and gluon GPDs in LO and NLO, 
using MRSTOl input PDFs, at the input scale Qo = 1 GeV 
for ( — 0.1 (upper plot) and Cj" = 0.001 (lower plot), values 
typical of HERA and HERMES kinematics, respectively. 



the input scale of CTEQ6M, as done below, we find an 
appropriate reduction in the cross section at low much 
more in line with the low data and the AJM value. 

It is important to note that the preliminary ZEUS data 
lies systematically above the HI data (see Fig. 11 of [^). 
Overall NLO seems to be doing better than LO, partic- 
ularly on the slope of the energy dependence. It is fair 
to say that all of the theory curves appear to have a Q^- 
dependence that is too steep to describe all of the data. 
We will return to this point in the next section. 

The difference between the MRST and CTEQ curves 
at LO and NLO reflects the relative size of the quark 
singlet and gluon distributions for each set. It is possible 
that more precise data on DVCS may eventually allow 
a discrimination between various input scenarios using 
NLO QCD. For this to be realistic one would first need 
to pin down the uncertainty associated with the slope by 
explicitly measuring the i-dependence. 

We also investigated the effect on the cross section 
of increasing the input scale for skewed evolution us- 
ing CTEQ input distributions, from the starting scale 
Qo = 1.3 GcV to Qo = 2.0 GeV. We then use the for- 
ward PDFs at the new scale in our model for the GPDs. 
Fig. |ll| shows that the reduced lever arm for skewed evo- 
lution starting at the higher scale leads to a smaller cross 
section at LO and NLO, as expected, and that, in LO 
at least, the effect of this change is rather large. In fact 
the CTEQ and MRST collaborations only advocate the 
use of their forward PDFs above w 3 - 4 GeV^ (they 
start evolution at a lower scale w 1 — 2 GeV^ due 
to technicalities associated with a consistent implemen- 
tation of charm). Hence, it is not completely clear where 
one should start skewed evolution, and this constitutes 
an additional uncertainty in the theoretical predictions. 

Having compared to small Xbj data we will now test 
the AJM ansatz for large Xbj by comparing to data on 
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FIG. 9: The photon level cross section, (j{"i*p IP), cal- 
culated using the forward model ansatz for input GPDs, in 
the average kinematics of the HI data: as a function of W at 
fixed — 4.5 GeV^ (upper plot), and as a function of at 
fixed = 75 GeV (lower plot). A constant slope parameter 
of -B = 6.5 GeV^'^ was used. 



the single-spin asymmetry, SSA, (HERMES and CLAS 
[|l9|) and the charge asymmetry, CA, (HERMES only), 
defined by 



SSA 



2/q d0 sm{(t)){da^ -da^) 
2 J^^ d(j) cos(0)((ifT+ — da^} 



J^^ d(j) {da~^ 



da 



(27) 



Here da\ da^ refer to the differential cross sections with 
the lepton polarized along or against its direction of mo- 
tion, respectively; da^,da^ are the unpolarized differ- 
ential cross sections for positrons and electrons, respec- 
tively. 

Such a comparison of QCD models with the avail- 
able high Xbj data may be viewed with some scepticism, 
especially in the case of CLAS which has such a low 
Q2 ^ 1 - 2 GeV^ (HERMES is only shghtly better with 
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FIG. 10: The photon level cross section, u{'y*p — » 7p), cal- 
culated using the forward model ansatz for input GPDs, in 
the average kinematics of the preliminary ZEUS data: as a 
function of W at fixed — 9.6 GeV^ (upper plot), and as a 
function of at fixed W = 89 GeV (lower plot). A constant 
slope parameter of _B = 6.5 GeV~^ was used. 



FIG. 11: The effect of changing the starting scale, Qo, on 
photon level cross section (t(7*p — > 7p), calculated using the 
forward model ansatz and CTEQ input PDFs, in the average 
kinematics of the preliminary ZEUS data: as a function of 
W at fixed = 9.6 GeV^ (upper plot), and as a function 
of at fixed W — 89 GeV (lower plot). A constant slope 
parameter of _B = 6.5 GeV~^ was used. 



a typical Q'^ of '--^ 2 — A GeV^). Firstly, it is a priori not 
clear that perturbation theory is applicable at such low 
values, (in particular, higher twist corrections may 
be expected to become important in this region and our 
approximations correspond to the DVCS cross section 
being divergent as 0.). Secondly, the previously 

neglected GPDs H, E and E become increasingly impor- 
tant as X increases [|l||§. In the following we will in- 
clude the dominant twist-3 contributions pdj , which are 
entirely kinematic in origin, in our calculation of the dif- 
ferential cross section, neglecting the sub-dominant twist- 
3 effects. We use the same input models for i/, E and E 
as well as i-dependence for the various amplitudes i.e. 
various dipole form factors which are about equivalent 
with an exponential at small t with a slope B between 5 
as in 1 21 



ourselves to MRSTOl input PDF for simplicity. 

It transpires that the average kinematics of HERMES 
is such that H is still the leading GPD and within our 
model assumptions i/, E and E could be set to zero for 
those values, with negligible difference to the final answer 
psf . Within the above caveats, we find for average HER- 
MES kinematics {{x) = 0.11, (Q^) = 2.56 GeV^ (t) = 
-0.265 GeV^) 

SSA = -0.28 (LO), - 0.23 (NLO) 
CA = 0.12 (LO), 0.09 (NLO) (28) 

compared to the quoted experimental results pil 



and 8, 



SSA= -0.21 ±0.08 
CA = 0.11 ±0.07. 



For HERMES we perform a full LO and 
NLO QCD analysis, whereas for CLAS we are restricting 
ourselves to LO, i.e., we are testing handbag dominance 

with no or little evolution. Furthermore, we shall restrict For the average CLAS kinematics {{x) — 0.19, ((J'^) 



(29) 
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1.31 GeV^ {t) = -0.19 GeV^) we find 

SSA = 0.2 (LO) (30) 
compared to the experimental value (second reference of 

SSA = 0.202 ±0.041. (31) 

This demonstrates that the AJM ansatz works surpris- 
ingly well even at large Xbj giving us confidence in the 
AJM-based model and suggesting that a fit to the avail- 
able data should be possible without tuning too many 
input parameters. 

IV. A SIMPLE MODEL FOR THE SLOPE 
PARAMETER 

It was pointed out in that the i-slope of the DVCS 
cross section at small x should depend strongly on in 
the transitional region from of a few GeV^ to large 
Q^. At ^ 2 GeV^ it is natural to expect that the 
slope will be pretty close to that for exclusive p-meson 
production: B ^ 8 GeV~^ |45). For large the dom- 
inant contribution is governed by evolution trajectories 
which, at the resolution Qq^2 GeV^, originate from the 
gluon field. Hence we expect that in this case the slope 
will be given by the square of the two-gluon form factor 
of the nucleon a,t X,X — ( ^ Xbj ■ Recently it was 
demonstrated that for Xhj > 0.05 this i-dependence can 
be approximated in a wide range of t as 1/(1 — t/m2g)* 
with m^g ~ 1.1 GeV^. This corresponds to a t-slope of 
B ^ 3 GeV^^ for exponential fits j4^. At smaller x an in- 
crease of the slope is expected which could originate from 
several effects, including Gribov diffusion. Hence for the 
highest point of ZEUS of about 90 GeV^ we expect 
B = 3.5 ±0.5 GeV"^ The recent HI and ZEUS p-meson 
production data, for a PU-range similar to the DVCS ex- 
periments, indicates that the slope of p-production both 
for (Ti and ctt drops rather rapidly with increasing 
reaching S - 5 at ~ iq QcV^ [g. 

A simple parameterization which reflects the discussed 
constraints for the range of 2 < < 100 GeV^ is 

BiQ') = Bo(l-Cln(^^^'^ (32) 

with Bo ^ 8 GeV-2, Qa ^ 2 GeV^, C = 0.15 be- 
ing reasonable values for the various parameters. This 
gives B(Q2 = 9.6) = 6.1 GeV'^ and B(Q2 = 4.5) = 
7.0 GeV-2 at the average values of the ZEUS and HI 
data, respectively (in broad agreement with our chosen 
constant value of i? = 6.5 GeV~^). Fig. ^illustrates the 
effect of introducing this simple model on the description 
of the Q^-dependence of the ZEUS data. 

This modification of the B slope gives a great improve- 
ment in comparison with the data and shows how impor- 
tant an experimental determination of the B slope is. 
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FIG. 12: The effect on the DVCS cross section, in the av- 
erage kinematics of the ZEUS data, of introducing our sim- 
ple Q^-dependent model of Eq.(^) for B, the slope of the 
t-dependence. 

since it constitutes a large theoretical uncertainty at this 
point. 



V. CONCLUSIONS 

Using sound phenomenological guidelines such as the 
aligned jet model, we have constructed a model for gen- 
eralized parton distributions at the input scale. Within 
certain theoretical uncertainties (such as the exact shape, 
the input scale and the functional form of i? in VF and 
Q^) this model can be used in a NLO QCD analysis to 
describe the recent DVCS data from HI, ZEUS, HER- 
MES and CLAS within their experimental errors. In con- 
structing this model we have given a simple and flexible 
algorithm which can be easily incorporated into a fitting 
procedure. 

We have also demonstrated and explained the failure of 
the most widely used model for generalized parton distri- 
butions, the factorized double distribution based model, 
to describe the available DVCS data, when rigorously ap- 
plied in a LO or NLO QCD analysis in its current form. 

The modeling of the input GPDs is now sufficiently 
advanced to justify attempting to fit some of the input 
parameters directly to the available data. A basic analy- 
sis of the data would seem to favor a ^-dependence with 
a slope parameter, B, that depends on Q^. Hence, an 
accurate measurement of this slope is of crucial impor- 
tance for further progress of the comparison of theory 
and experiment. 
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